Mon. Not. R. Astron. Soc. 000,[T]l5](2009) Printed 22 September 2009 (MN I^TbX style file v2.2) 



Deciphering solar turbulence from sunspots records 



F. Plunian\ G. R. Sarson^, R. Stepanov^ 

' LGIT, UJF. CNRS. B.P. 53. 38041 Grenoble Cedex 9, France 

' School of Mathematics and Statistics, Newcastle University, Newcastle NEl 7RU, UK 
^ Institute of Continuous Media Mechanics, Korolyov 1, 614013 Perm, Russia 



Accepted 2009 .... Received 2009 in original form 2009 



ABSTRACT 

It is generally believed that sunspots are the emergent part of magnetic flux tubes in the solar 
interior These tubes are created at the base of the convection zone and rise to the surface 
due to their magnetic buoyancy. The motion of plasma in the convection zone being highly 
turbulent, the surface manifestation of sunspots may retain the signature of this turbulence, 
including its intermittency. From direct observations of sunspots, and indirect observations of 
the concentration of cosmogenic isotopes '''C in tree rings or "'Be in polar ice, power spectral 
densities in frequency are plotted. Two diff'erent frequency scalings emerge, depending on 
whether the Sun is quiescent or active. From direct observations we can also calculate scaling 
exponents. These testify to a strong intermittency, comparable with that observed in the solar 
wind. 
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I INTRODUCTION 

Sunspots observed at the surface of the convection zone of the Sun 
are usually understood as the manifestation of solar magnetic ac- 
tivity. Naked-eye and telescope observations of sunspots are avail- 
able from AD 1610, providing reliable records of sunspot num- 
bers (SSN). Several sets of data exist, varying in how they have 
been sampled, averaged (daily or monthly), whether they concern 
sunspots or sunspots groups, and on the scientific societies who 
have conipiled the records. Here we consider the American daily 
SSN (D)Q, the American daily group SSN (G)B and the Interna- 
tional monthly averaged SSN (M)|3 (figure[Tll. In addition SSN at 
earlier times have been reconstructed from proxies, based on the 
concentration of cosmogenic isotopes in tree rings or '"Be 
in ice core bubbles. The production rate of such isotopes increases 
with the cosmic-ray flux, which is higher when the solar magnetic 
activity is low. Plotting the SSN versus time reveals a cycle of about 

II years known as the Schwabe cycle. This cyclical solar mag- 
netic activity is sufficiently robust to be detected in '"Be concen- 
tration records, even during some long periods with al most no vis- 
ible s unspots, like the Maunder minimum (1645-1715) l lBeer et al] 
Il998h . Analyzing long time series of '"•€ and '"Be, it has been 
shown that the solar activity of the last 70 years has been excep- 



' ftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/SUNSPOTJ4UMBERS/ 
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tionally high jUsoskin et aLll2003l : ISoTanki et alj|2004l) . and that a 
decline is expected within the next two or three cycles l lAbreu et al.l 
120081) . 



The occurrence of sunspots is of course an important diag- 
nostic that must be reproduced by any solar dynamo model. The 
fact that it is irregular (in spite of the Schwabe cycle) reflects the 
complexity inherent in the nonlinear coupling between the tur- 
bulent flow and the m agnetic field in the solar convection zone 
terowning et alj2006l) . With the long time series of direct SSN ob- 
servations and '*C and '"Be data available, it is tempting to calcu- 
late the corresponding frequency spectra, to infer some signature 
of the underlying turbulence in the convection zone. Similar at- 
tempts have been made for the Earth jCourtillot & Le Moueil 19881 : 
[Constable & lohnsonl I2OO5I : ISakuraba & Hamanol |2007|) : the fre- 
quency spectrum of the geomagnetic dipole moment obtained from 
paleomagnetic data is consistent with an inertial range of /"^'' scal- 
ing, and a dissipatio n range of /~"^^ as ex pected in magnetohydro- 
dynamic turbulence l lAlemanv et al.l2000l) . For the Su n, analysis of 
International daily SSN hav e led to a /"^'^ scaling jMorfill et al.l 
ll99lULawrence et alj[l995l) , and this has been attributed to some 
sequential sampling of the field upon arrival at the photosphere on 
top of a Kolmogorov spatial scaling due to the underlying turbu- 
lence. Here we extend this analysis to the other SSN records men- 
tioned above and test how robust the /"^'^ scaling is, in particu- 
lar during minima and maxima of sunspot activity. In addition, for 
time scales smaller than 2 years, the stoc hastic character of the SSN 
records suggests strong intermittency jLawrence et al.l [l995h . as 
opposed to a low-dimensional chaotic interpretation. We shall char- 
acterize this intermittency by calculating the corresponding scaling 
exponents. 
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Figure 1. Data sets of 10 year-averaged SSN from '''C (top blue) and '"Be 
(top red), International SSN monthly averaged M (second), American daily 
group SSN G (third), and American daily SSN D (bottom). 
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Figure 2. Wavelet coefficients (logarithm of absolute value) of '*C (top), M 
(second), G (thkd), and D (bottom). 



2 WAVELET SPECTRA 

In order to filter out the noise, we use a wavelet decomposition of 
the signal. In figure (2] the wavelet coefficients are plotted versus 
time (in years) and frequency (in month"'). The light (resp. dark) 
colors coffespond to low (resp. high) values of these coefficients. 
For a given year, the curve giving the wavelet coefficient versus 
frequency corresponds to a power spectral density of the signal. 
The dark horizontal stripe for / ~ 0.01 month"' corresponds to 



the Schwabe cycle. It is not visible in the C data set due to the 
coarse sampling of this data (averaged over 10 years). On the other 
hand a clear dark stripe is visible for / ~ 3.5 x 10"^ month"', cor- 
responding to a ~2400 years cycle. In the figure for M data the 
dark stripe of the Schwabe cycle almost disap pears during the Dal- 
ton minimum (1790-1820) l lFrick et al.ll 19971) . In the G and D data 
figures we identify another horizontal stripe at / ~ 1 month"' , cor- 
responding to the solar mean rotation rate. This indicates that the 
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Figure 3. Wavelet spectral density versus frequency, on average (top), and 
for magnetic activity minima (bottom, dashed curves) and maxima (bottom, 
solid curves). The dashed straight hnes correspond to —2/3 and -8/3 slopes 
(top), and -2/3 and —1 slopes (bottom). 



latitudinal repartition of sunspots is not homogeneous. Finally the 
vertical light stripes correspond to minima of magnetic activity. 

The time average of the wavelet coefficients are shown in fig- 
ure[3](top) for the five data sets. As mentioned earlier the two peaks 
around 0.01 month"' and 1 month"' correspond to the Schwabe cy- 
cle and solar mean rotation rate. Between them the three spectra of 
the directly observed SSN (M, G, D) present a common scaling in 
/~^'^. The other data sets '''C and "'Be are compatible with a /""'^ 
scaling as well, although the ''*C PSD is overestimated by roughly 
a factor 10, probably due to the proxy used in the reconstruction of 
the SSN from the '''C data. 

For the three data sets (M), (G) and (D), instead of averag- 
ing on all times, we now average on periods corresponding to ei- 
ther maximum or minimum magnetic activity as shown in figure 
|4l For the maximum (resp. minimum) activity subset, the excluded 
data is that centred around the times of the Schwabe minima (resp. 
maxima) The corresponding spectral densities are plotted in fig- 
ure [3] (bottom). The slopes for minima are systematically steeper 
than those for maxima, indicating two different regimes. To esti- 
mate these slopes we vary both the range of frequency /inax] 
on which they are calculated, and the way the data sets are split 
into subsets of maximum and minimum activity. For the former 
we take = 0.7 month"' for the data sets (G) and (D) in or- 
der to escape from the influence of the peak / = 1 month"', and 
/max = 0.5 month"' for the data set (M). When changing the 
slopes change. We vary such that the ratio fmax/fmm is about 
10, and the standard deviation of the slope remains 10 % or less 
of its average value. This leads to /„i„ € [0.05,0.07] for (G), 
/m,„ € [0.06, 0.08] for (D), and /^,„ € [0.03, 0.05] for (M). In addi- 
tion we consider at least three different degrees of splitting for each 
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Figure 4. Same data as in figure[T](the last three panels) but subdivided into 
sets of maximum (dark) and minimum (light) magnetic activity. 



data set, this splitting degree being related to the time length of 
the subsets of maximum and minimum activity. The choices of this 
splitting degree are made such that both subsets are long enough 
to provide good statistics, but remain separated by sufficient time 
lags that the SSN values falling into each subset do not overlap 
too much. Then for both the maximum and minimum subsets, we 
calculate the mean slope and the standard deviation obtained when 
varying both the frequency range and the degree of splitting. The 
corresponding slope estimates are given in table[T] They are consis- 
tent with power spectra in /"^'^ and /"' . The corresponding dashed 
lines are plotted in figure[3]to guide the eye. The standard deviations 
are small, showing that these slopes are robust with respect to the 
details of our analysis. The formal standard errors from each of the 
individual regressions (for specific frequency ranges and degrees of 
splitting) are of comparable magnitude. 

As noted by Lawrence et al. (1995), the question of causality 
complicates the interpretation of such temporal data. The difference 
of spectral slopes between minima and maxima can be attributed 
to two different effects: a change of the underlying turbulence, af- 
fecting the spatial structure the magnetic field; or a change in the 
frequency of the sequential sampling of the magnetic field, as sug- 
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Activity 

Max. 
Min. 



(M) 



(G) 



-0.61 ±0.05 -0.69 ±0.05 
-1.03 ±0.03 -0.94 ±0.07 



(D) 

-0.63 ± 0.05 
-0.98 ± 0.02 



Table 1. Slope estimates for the PSD curves plotted in figure[3](bottom). 
They correspond to average values plus standard deviation errors when 
varying both the frequency range and the degree of splitting of each data 
set (into the two subsets of maximum and minimum magnetic activity). 



gested by Lawrence et al. (1995). Although the latter effect cannot 
be excluded, there is a simple argument in favour of the former. It 
is generally accepted that the occurrence of sunspots at the photo- 
sphere is due to the magnetic buoyancy force V6', wh ere B is some 
magn etic induction intensity in the convection zone jTobias et al.l 
I2OOII) . It is then tempting to interpret the two spectral slopes as the 
signatures of this buoyancy, assuming that the frequency of sunspot 
occurrence at the photosphere is proportional to this force. Then 
the /"^'^ and SSN spectra would correspond to buoyancy spec- 
tra of kr'^^^ and kr^ , where k is the spatial wave number. During 
maxima this implies a Kolmogorov magnetic energy spectrum of 
compatible with inertia-driven turbulence in the convection 
zone. During minima it implies a magnetic energy spectrum of k^^, 
comp atible with turbulence dominated by the solar rotation jZhoul 
1 19951) . In the transport scenario proposed by Tobias et al. (2001), 
the field which arises at the surface is the strongest part of a poloidal 
field generated by cyclonic turbulence in the convection zone. Our 
interpretation then suggests two different regimes for such cyclonic 
turbulence, controlled by either inertia or rotation. 



3 INTERMITTENCY 

The stochastic nature of SSN occurrence for times scales smaller 
than 2 years has been shown by Lawrence et al. (1995), suggesting 
an intermittent turbulence. Here our goal is to quantify this inter- 
mittency for the three sets (M, G, D), calculating the corresponding 
scaling exponents. For that we first calculate the associated gener- 
alized structure function (GSF) and look for its scaling exponents, 
as usually done in turbulence. We define the SSN increment by 



6y{t,r) = S(t + T)-S{t) ., 



(1) 



where 5 (f) denotes the SSN at time t. Assuming statistical station- 
arity in the frequency range of interest, the t dependence in 6y(t , t) 
can be dropped and the GSF is then given by jNicol et alJ2008h 



s,„(T) = m'"} 



\6y\"' P{Sy,T)d(Sy) 



(2) 



where P is the probability density function (PDF) of 6y, the angle 
brackets (■) denote time averaging, and m is a positive integer 

In figure |5] PDFs for the three data sets (M, G, D) are plot- 
ted for selected values of t. The PDFs of the '^C data are poorly 
defined, and so we drop this data set for the rest of the study. The 
PDFs of the other data sets show peaks at t = 11/2 years, corre- 
sponding to the Schwabe cycle. For other values of t they exhibit 
tails containing a higher number of rare events than for a gaussian 
distribution, suggesting intermittency. Similar results were shown 
in Lawrence et al. (1995). 

To quantify this intermittency we first check whether the GSF 
obey a scaling law in the form 



SjT) 



(3) 



We find (not shown) that this is clearly the case. In homogeneous 
and isotropic fully developed turbulence, intermittency corresponds 
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Figure 5. Probability density functions oi S{t + t) - S{t) for M (top), G 
(middle) and D (bottom). The labels indicate the value of r in years. 



to ({m) < m/3. The ratio ((m)/((3) is calculated, estimating the 
scaling power of 5„,(T)/5m(3). Plotting the ratio ({m)/^(3) versus 
m for the three data sets (figure[6]l we see a clear departure from the 
Kolmogorov straight line ^(m)/((3) = m/3, and clear indications 
of intermittency. It is remarkable that the (D) set, which has the 
best sampling, leads to the largest intermittency. It is also remark- 
able that the (M) and (G) sets lead to similar scaling exponents, 
supporting the equivalence between averaging over space and time. 

The exponent s can be fitted to the standard p -model derived 
for hydrodynamic jMeneveau & Sreenivasaniri987h and magneto- 
hydrodynamic turbulence l lCarbonelll993h 7This model is defined 
by 



<r(m)= l-log2[p'"/3+(l-p)'"/3j. 



(4) 



We find P(G) = 0.68, P(M) = 0.68, and = 0.83. The last value 
compares surprising wel l with those for the solar wind measured b y 
the Ulysses spacecraft jPagel & Balogi]|2002l : InIcoI et al.ll2008l) . 
and for the magnetosp heric cusp measured by the Polar satellite 
jYordanova et al.l2004h . even though the frequencies differ by sev- 
eral orders of magnitude. 
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Figure 6. Scaling exponents l^(m)/({3) plotted versus m for the three sets 
of data M, G, D, and for 2 months < r < 14 months. The dashed line 
corresponds to a Kolmogorov scaling f (m)/f (3) = m/3. 
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4 SUMMARY 

In conclusion, the wavelet spectral analysis of sunspot records has 
revealed two different behaviors, depending on whether the Sun is 
quiescent or active. This suggests two different kinds of turbulence 
in the convection zone, controlled either by inertia or by rotation. 
The signature of such fully developed turbulence is confirmed by 
the calculation of the GSF scaling exponents, which indicate strong 
intermittency. 
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